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Abstract 



In the paper, the large time behavior of solutions of the Cauchy problem for the 
one dimensional fractal Burgers equation u t + {—d 2 ) a / 2 u + uu x = with a £ (1, 2) 
is studied. It is shown that if the nondecreasing initial datum approaches the 
constant states u± (u_ < u + ) as x — ► ±oo, respectively, then the corresponding 
solution converges toward the rarefaction wave, i.e. the unique entropy solution of 
the Riemann problem for the nonviscous Burgers equation. 
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1 Introduction 

The goal of this work is to study asymptotic properties of solutions to the Cauchy problem 
for the nonlocal conservation law 



where A" = (—d 2 /dx 2 ) a ^ 2 is the pseudo differential operator denned via the Fourier trans- 
form 



u t + A' 



u + uu x = 0, x e R, t > 0, 



(1.1) 





(A<*v)(0 = \t\ a v(0- 
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Following [2J, we will call (jl.ip the fractal Burgers equation. Equations of this type 
appear in the study of growing interfaces in the presence of selfsimilar hopping surface 
diffusion |17j . Moreover, in their recent papers, Jourdain, Meleard, and Woyczynski 
[EH EE2] gave probabilistic motivations to study equations with the anomalous diffusion, 
when Laplacian (the generator of the Wiener process) is replaced by a more general 
pseudodifferential operator generating the Levy process. In particular, the authors of [12] 
studied problem (ll.lj) - (ll.2p . where the initial condition uq is assumed to be a nonconstant 
function with bounded variation on R. In other words, a.e. on R, 

/X 
m(dy) = c + H * m(x) (1.4) 
-oo 

with c G R, m being a finite signed measure on R, and H(y) denoting the unit step 
function 1{<,>o}- Observe that the gradient v(x, t) = u x (x,t) satisfies 

v t + A a v + (vH * v) x = 0, u(.,0) = m. (1.5) 

If m is a probability measure on R, the equation (11.51) is a nonlinear Fokker-Planck 
equation. In the case of an arbitrary finite signed measure, the authors of (12] associate 
(11. 5p with a suitable nonlinear martingale problem. Next, they study the convergence of 
systems of particles with jumps as the number of particles tends to +oo. As a consequence, 
the weighted empirical cumulative distribution functions of the particles converge to the 
solution of the martingale problem connected to (11.51) . This phenomena is called the 
propagation of chaos for problem (jl.ip - fll.2D and we refer the reader to [12] for more 
details and additional references. 

Motivated by the results from [12], we study problem (jl.ip - fll.2l) under the crucial 
assumption a G (1,2) and with the initial condition of the form (jl.4p . In our main result, 
we assume that uq is a function satisfying 

Uq — u_ G L ((— oo, 0)) and Uq — u + G L ((0, +oo)) with n_ < u + , (1.6) 

where w_ = c and u + = c + f R m(dx). 

It is well known (cf. [TUj [TSl [H] and Lemma 12.41 below) that the asymptotic profile as 
t — > oo of solutions to the viscous Burgers equation 

u t -u xx + uu x = (1.7) 

(i.e. equation (jl.lj) with a = 2) supplemented with an initial datum satisfying (j 1.6 ft is 
given by the so-called rarefaction wave. This is the continuous function 

i- , x/t <U-, 

w R (x,t) = W R (x/t) = { x/t, u„<x/t<u + , (1.8) 

u + , x/t > u+, 

which is the unique entropy solution of following Riemann problem 

w R + w R w R = 0, (1.9) 
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w R (x, 0) = w§(x) 



{ 



U- , x < 0, 
u + , x > 0. 



(1.10) 



Below, we use the solution of the Burgers equation (11.71) with the initial datum (11.101) as 
the smooth approximation of the rarefaction wave (11.81) . 

The authors of this work were inspired by the fundamental paper of Il'in and Oleinik 
[TO] who showed the convergence toward rarefaction waves of solutions to the nonlinear 
equation u t — u xx + f(u) x = under strict convexity assumption imposed on /. That 
idea was next extended in several different directions and we refer the reader e.g. to 
[HI [HI [19j [2U [22] for an overview of know results and additional references. 

In this work, we contribute to the existing theory by developing tools which allows 
to obtain analogous results for equations with a nonlocal and anomalous diffusion. Basic 
properties of solutions (namely, their existence and the regularity) of quasilinear evolu- 
tion equations with (— A) Q / 2 , a G (1, 2), (or, more generally, with the Levy diffusion) were 
shown in [TIE]- On the other hand, one may expect singularities in finite time of solutions 
to (11.11) with a G (0, 1), see [TJ for more details. If uo G i/iTSL), Biler, Karch, and Woy- 
czynski [3] proved that the large time asymptotics of solutions to (jl.lj) - (11.2D is described 
by the self-similar fundamental solution of equation v t + A a v = 0. Analogous asymptotic 
properties of solutions to multidimensional generalizations of problem (jl.lD - fll.2p with 
u G L\R N ) are studied in [HE]. 

Finally, we would like to report the recent progress in the understanding of properties 
of solutions of the quasi-geostrophic equation with an anomalous diffusion, cf. [HI [13] 
and the references therein. We are convinced that our techniques can be applied to that 
equation, as well. 

The purpose of the present paper is to prove the convergence of solutions to the Cauchy 
problem for the fractal Burgers equation (jl.ip - fll.2fl toward rarefaction waves. We state 
our main result in the following theorem. 

Theorem 1.1. Let a G (1,2). Assume that w R = w R (x,t) is the rarefaction wave (I l.S\) . 
Denote by u = u(x,t) the unique solution of problem (jj. jp - flj.jjp corresponding to the 
initial datum Uo of the form l \1.4\) and satisfying (j 1.6§ . For every p G ((3 — a)/ (a — 1), oo] 
there exists C > independent oft such that 



for all t > 0. 

Remark 1.1. Our result and its proof hold true also for a = 2 (observe that (3 — 
a) /{a — 1) — > 1 as a — ► 2). However, we pass over this case for simplicity of the exposition 
and because the large time asymptotics of solutions to the Burgers equation (11.71) is well- 



In the next section, we gather several preliminary properties of the operator A a and 
of solutions to problem (jl.lfl - fll.2p . Theorem 11.11 is shown in Section 3. In Section 4, we 
discuss possible generalizations of our main result. 



u(t) — w 



R 



(t)\\ P < cr^- 1 ^ 3 -^ 2 io g (2 + 1) 




□ 
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Notation. For 1 < p < oo, the L p -norm of a Lebesgue measurable, real-valued 
function v defined on R is denoted by ||u|| p . For a finite signed measure m on R, we 
put ||m|| = |to|(R), where \m\ is the total variation of m. The Fourier transform of v 
is v(£) = (27f)^ 1 ^ 2 J R e~ lx ^v(x) dx. Given a function v = v(x), we are going to use the 
decomposition v — v~ + v + , where as usual v~ = max{0, —v} and v + = max{0, v }. The 
constants (always independent of of t) will be denoted by the same letter C, even if they 
may vary from line to line. Occasionally, we write, e.g., C = C(a,£) when we want to 
emphasize the dependence of C on parameters a and £. 



2 Preliminary results 

We begin by recalling that the basic questions on the existence and the uniqueness of 
solutions of problem (jl.ip - fll.2D were answered in the papers OH]. 

Theorem 2.1. ([7, Thm. 1.1], jgj Thm. 7]) Let a G (1, 2) and u G L°°(R). There 
exists the unique solution u = u(x,t) to problem fl 1. iD - fl Lift) in the following sense: for all 
T > 0, 

u G C 6 ((0, T) x R) and, for all a G (0, T), u G C 6 °°((a, T) x R), 

u satisfies dOD on (0, T) x R, 

u(t, •) — >■ m in L°°(R) weak — *, as t — > 0. 

Moreover, the following inequality holds true 

IK*)IU < IKHoo for all t>0. (2.1) 

The main goal of the section is to complete this result by additional properties of u x 
if the initial conditions are of the form (11.41) . 

Theorem 2.2. Let a G (1, 2). Assume that the initial datum uq can be written in the 
form fll.^p for a constant c G R and a signed finite measure m on R. Then the solution 
u = u(x,t) of problem (I T, Oj) satisfies u x G C((0, T]; L P (R)) /or eac/i 1 < p < cxd and 
every T > 0. 

Consider u and u two such solutions with initial conditions u and u , respectively. 
Suppose that u x (x,t) is nonnegative a.e. and u — u G L X (R). Then 

||u(0-m(0I|i< K-Wolli (2-2) 

/or a// 1 > 0. 

Theorem 2.3. Under the assumption of Theorem \2.^ if the measure m in the initial 
datum l \1.4\ ) is nonnegative, we have 



i) u x (x, t) > for all x G R and t > 0, 

/or every p G [1, oo] there exists C = C(p) > such that 



< r 1+1/p \\m\\ 1/p . (2.3) 
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In the proofs of Theorems 12.21 and 12.31 as well as in our study of the large time asymp- 
totics, we shall require several properties of the operator A a and of the semigroup of linear 
operators generated by it. First of all, note that the operator defined by (11.31) has the 
integral representation for every a £ (1, 2) (cf. eg. [8j Thm. 1]) 

A°w(z) = -C(a) [ ^ + dz {2A) 

Jr \ z \ 

This formula allows us to apply A a to functions which are bounded and sufficiently smooth, 
however, not necessary decaying at infinity. 

Lemma 2.1. Let 1 < a < 2. For every p £ [1, oo] there exists C = C(p,a) > such that 

IIA^II^CII^II^II^H-- 1 (2.5) 
all functions w satisfying w x ,w xx £ L P (IR). 

Proof. We can easily deduce the interpolation inequality (I2.5P from (|2.4p . Indeed, it 
follows from the Taylor formula that for any fixed R > we have 

||A a u>|| p < Cllm^Hp / \ z \ 1 ~~ a dz + C||wa;||p / \z\~ a dz 

J\z\<R J\z\>R 

< C(R 2 - a \\w xx \\ p + R l - a \\w x \\ p ). 

Choosing R = \\w x \\ p /\\w xx \\ p we complete the proof of inequality (12.51) . □ 

Now, we prove the Nash inequality for the operator A a . 
Lemma 2.2. Let < a. There exists a constant > such that 

\H\f +a) <C N \\A a / 2 w\\l\H\l a (2-6) 
for all functions w satisfying w £ L 1 (M) and A a ^ 2 w £ L 2 (R). 

Proof. For every R > 0, we decompose the L 2 -norm of the Fourier transform of w as 
follows 



\H\ 2 2 = c / 

Jr 

<CNI«/ dt + CR- a [ \£\ a \w(0\ 2 dZ 
J\Z\<R J\£\>R 

< CR\\w\\ 2 + CR~ a \\A a/2 w\\j. 
For R = (||A Q / 2 ^||2/||^|| 2 ) 1/(1+Q) we obtain (Q. □ 
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{A a w)\w\ p ~ 2 w dx > ^ P 2 ^ / (^\w\i) 2 dx (2.7) 



Lemma 2.3. Let < a < 2. For every p > 1, we have 

4(p- 1) / ^ A a,_ ,g\2 

for all w G L P (R) suc/t i/m* A Q w G L P (M). // A Q w G ^(R), we obtain 

{A a w)sgnwdx>0, (2.8) 

K 

and, ifw,A a w G L 2 (R) ; it follows 

I {A a w)w + dx>0 and [ {A a w)w~ dx > 0, (2.9) 
Jr Jr 

where w + = max{0, w} and w~ = max{0, — w}. 

Inequality (12. 7p is well-known in the theory of sub-Markovian operators and its state- 
ment and the proof is given e.g. in [T5l Theorem 2.1 combined with the Beurling-Deny 
condition (1.7)], see also [6] [13]. Observe that if a = 2, integrating by parts we obtain 
(12. 7ft with the equality. Inequality ( 12.81) (called the Kato inequlity) is used in [8] to con- 
struct entropy solutions of (11. lft and it can be easily deduced from [8] Lemma 1] by an 
approximation argument (see also [2] Inequality (3.5)]). The proof of (12.91) can be found, 
for example, in [151 Proposition 1.6]. 

We also recall that, by Duhamel's principle, the solution to problem (ll.ll) - (ll.2p can be 
written in the equivalent integral form 

u(t) = S a (t)uo — / S a (t — t)u(t)u x (t) d tau, (2-10) 
Jo 

where 

S a (t)u =p a (t) * u (x). (2.11) 

Here, the fundamental solution p a (x,t) of the linear equation d t v + A a v = can be 
computed via the Fourier transform p a (£,t) = e~'^ Q . Hence, p a (x,t) = t _1 / a P a (xt _1 / a ), 
where P a is the inverse Fourier transform of e~'^' Q . It is well known that for every a G (0, 2] 
the function P a has the property L P a {x) dx = 1 and is smooth, nonnegative, and satisfies 

< P a {x) < C(l + \x\y {a+1) and \d x P a (x)\ < C{1 + |x|)~ (a+2) (2.12) 

for a constant C and all Using these properties of the convolution operator S a (t) 

defined by (12.1 ip we obtain the estimates 

\\S a (t)v\\ p < Ct-t-WMWvW!, (2.13) 
\\(S a (t)v) x \\ p < Cr( 1 - 1 ^/ a - 1 /°||t;||i (2.14) 

for every p G [1, oo] and all t > 0. Moreover, we can replace v in ( 12 . 131) and in (12. 14ft by 
any signed measure m. In that case, ||i>||i should be replaced by ||m||. 
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Proof of Theorem \2.2\ It follows from the integral equation (I2.10p that u x is the solution 
of 

u x (t) = S a {t)m - I d x S a (t - t)V(t)u x (t) air, (2.15) 
Jo 

where V(x,t) = u(x,t) is treated as given, smooth, and bounded. Now the standard 
argument involving the Banach fixed point theorem allows us to show that the solution of 
the "linear" equation (12.151) has the solution in C((0, T]; L P {R)) for each p G [1, oo] and 
every T > 0. Here, we should use the following estimate of the operator T(u) defined by 
the right-hand side of (12.151) 

\\T(u)(t)\\ p <\\S a (t)m\\ p + [ \\d x S a (t-T)V(T)u x {r)\\ p dr 

t 



o 



KCt-WMoWmW+C sup Hn^lU / {t-r)- l / a \\u x {T)\\ p dT 

t€[o,t\ Jo 



being the immediate consequence of (12.131) and (12.141) . Let us skip other details of this 
well-known argument (cf. [20]). 

Now, we prove inequality (12.21) . A direct calculation shows that the function v (x, t) = 

u(x, t) — u(x, t) satisfies 

v t + A a v + ^{v 2 + 2vu)) x = 0. (2.16) 
First, we multiply equation (12.161) by sgnt> = t>|t>| _1 : 

— [ \v\dx+ [ (A a t>)sgnt> dx H — j [v 2 + 2vu)] x sgnv dx = 0. 
dt J-r 2 J R 

The second term is nonnegative by (12. 8ft . To show the same property for the third term, 
we replace the sgn function by smooth and nondecreasing (p = tp(x). In this case, we 
obtain 

/ [v 2 + 2vu\ x (f(v)dx = — / (v 2 + 2vu)(p'(v)v x dx — — ^(v ) x dx + u x Q(v) dx, 
Jr Jr Jr Jr 

where ^f(s) = J Q S z 2 ip'(z) dz and $(s) = J Q S 2z<p'(z) dz. Obviously, the first term on the 
right hand side is equal to zero and the second one is nonnegative because u x > and 
$(s) > for all s G R Now, the standard approximation argument gives J R [v 2 + 
2vu\ x sgnv dx > 0. Hence ||f(t)||i = \\u(t) — u(t)\\i < \\u — u Q \\i = \\v \\i for all t > 0. □ 

Proof of Theorem \2.3[ To show part (i) of Theorem 12. 3[ we deal first with the smooth 
initial datum u Q satisfying u 0tX (x) > and u 0)X G L P {R) for every p G [1, oo]. In this case, 
differentiating equation (11.11) with respect to x we have 

(u x ) t + A a u x + (uu x ) x = 0. (2.17) 
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Note the well known property 



I'll' tlx = j t'l r tlx — — / (v ) 2 cfa\ 



v<0 * U,L Jv<0 



Hence, multiplying (12.171) by u x , integrating the resulting equation over IR, and integrating 
by parts on the right hand side, we obtain 



\t+ I M 2 dx+ [ (A a u x )u x dx 

A al Ju x <0 Jr 



[uu x ) x u x dx 

u x <0 



{u x fdx. 

u x <0 



Since J R (A a u x )u~ dx > by (12. 9p and f u <Q (u x (x, 0)) 2 dx = by the assumption imposed 
on m , the Gronwall inequality implies J u <Q (u~(x, t)) 2 dx = for all t > 0. Consequently, 
u x (x, t) = and the proof of (i) for regular initial conditions is finished. 

Now, the proof of (i) for the solution u = u(x,t) corresponding to the initial datum 
Mo of the form (11.41) with the nonnegative finite measure m can be completed by the 
following approximation argument. We consider the sequence of regular initial conditions 
Uq as in the first part of this proof. Moreover, we assume that u'q x converges weakly 
to m and \\u'q — u ||i — > as n — > oo. Inequality (12.21) allows us to prove that the 
corresponding solutions u n (-,t) satisfy \\u n (-,t) — u(-,t)||i — > as n — >• oo for any t > 0. 
Hence, there is a subsequence n k — >• oo such that u nk (x,t) — > u(x,t) a.e. Since each 
u n (x,t) is nondecreasing as function of x, the same conclusion holds true for u(x,t). 

In order to show inequality (12. 3p . we first observe that integrating equation (I2.15P over 
IR and using the equalities 

/ S a (t)m dt = m(dx) and / d x S a {t — t)(u(t)u x (t)) dx = 
Jr Jr Jr 

we obtain the identity j R u x (x,t) dx = J R m(dx) which for nonnegative u x means 

|K:(t)||i = ||m|| for all t > 0. (2.18) 

Now, for fixed p G (1, oo) and u x > 0, we multiply (12.171) by u v ~ x and integrate the 
resulting equation over R. After some manipulations involving integrations by parts on 
the right hand side we arrive at 



1 d | 
p dt 



u x\\p+ / u p x 1 A a u x dx = — / (uu x ) x u x 1 dx 
Jr Jr 



p — 1 
p 



[ u p x +1 dx. (2.19) 
Jr 



Recall now that J r m^ 1 A a u x dx > by inequality ( 12. 7ft . Moreover, it follows from the 
Holder inequality combined with ( 12.18P that 



I^WIlf^^KWIlKilHr^- 
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Applying those two inequalities to (12.191) (note u x > 0) we obtain the following differential 
inequality for ||iia;(0llp 

^iiu.(t)ig < -(p- ^iimir 1 /^ 1 ) (ii^wig)^^. 

Integrating it we complete the proof of (I2.3P for any p G (1, oo). 

The case of p = oo is obtained immediately passing to the limit p — > oo in inequal- 
ity Q. □ 

We conclude this section by recalling some results on smooth approximations of rar- 
efaction waves, namely, the solutions of the following Cauchy problem 

w t - w xx + ww x = 0, (2.20) 

{it— x <c 
(2.21) 
u+ , x > 0. 

Lemma 2.4. Let < u + . Problem (\2.2U\) - (\2.21\) has the unique, smooth, global-in-time 
solution w(x,t) satisfying 

i) u_ < w(t, x) < u + and w x (t, x) > for all (x, t)6lx (0, oo); 

ii) for every p G [1, oo], there exists a constant C = C(p,u^, u + ) > such that 

\\w x (t)\\ p < cr 1+1 / p , \\w xx (t)\\ p < ct-v 2+l '^ 

and 

\\w(t) -w R (t)\\ p < Ct-^ 1M/2 , 
for all t > 0, where w R (x,t) is the rarefaction wave (I l.S\) . 

All results stated in Lemma 12.41 are deduced from the explicit formula for solutions to 
(I2.20p - (l2.2ip and detailed calculations can be found in [9] with some additional improve- 
ments contained in [TU Section 3]. 

3 Convergence toward rarefaction waves 

For simplicity of the exposition, we divide the proof of Theorem 11.11 into a sequence of 
Lemmata. 

Lemma 3.1. Let a G (1,2). Assume thatu andu are two solutions of problem 

with initial conditions uq and uq, the both of the from (| 1.$ with finite signed measures 



m and m, respectively. Suppose, moreover, that the measure rh ofu is nonnegative and 
uq — u G L 1 (R). Then, for every p G [1, oo] there exists a constant C = C(p) > such 
that 

\\u(t) - u(t)\\ p < Ct- {1 - 1/p)/a \\u - u \\i (3-1) 

for all t > 0. 
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Proof. In our reasoning, we denote v(x,t) = u(x,t) — u(x,t) which satisfies equation 
(gJBD - It follows from Theorem inequality (Q, that ||f(t)||i < ||vo||i- 
Now, we multiply equation (I2.16P by |w| p_2 t' with p > 1 

™ / M P dx + f (A a v)(\v\ p - 2 v)dx + \ [ [v 2 + 2v^ x \v\ p ~ 2 v dx = 0. (3.2) 
P dt 7r 7r 2 7k 

The third term on the left hand side of ( 13. 2j) is nonnegative by the following calculations 

[v 2 + 2vu\ x \v\ p - 2 v dx = / 2^|t;| p rfx + / 2u^M p - 2 t; + / 2E c H P dx (3.3) 
</r 7r </r Jr 

= 211 1 / uJv\ p dx > 0, 



because J R v x \v\ p dx = and u x > 0. Hence, using inequality (12. 7p . we obtain from (13. 2p 



4- [ \v\ p dx + 4(l--] [ (A a/2 \v\ p/2 ) 2 dx < 0. (3.4) 
at 7r V P/ 7r 



From now on, we proceed by induction. Applying the Nash inequality (12.61) combined 
with ( 12.21) . we deduce from ( 13.41) with p = 2 the following differential inequality 



d 



^)ll 2 + 2^ 1 i| Vo |ir 2Q ll^)ll2 (1+a) <o, 



which, after integration, leads to 

||«(*)||2<Ci||i;o||ir 1/(2a) with d = (C N /2a) 1,{2a) . (3.5) 

This is estimate ( 13.11) with p = 2. 
Suppose now that 

\\v(t)\\ 2 n <Cn* _(1-2 "" )/a |ko||i for all t>0. (3.6) 

We consider (13.41) with p = 2 n+1 , where the second term is estimated, first, by the Nash 
inequality (12.61) with w = \v | 2 ™, next, by the inductive hypothesis (13. 6p . This two-step 
estimate leads to the differential inequality 

f t \wt)\\t: + + i +4(1 - 2- n - i )c N \c n \\v \\ 1 )- a2n+1 t 2n+l - 2 (\Htn 2 2^) 1+a < o. 

Integrating it we obtain 

IKt)|| 2 n+i < C n+1 r (1 - 2 ^ 1)/a ||i;o||i for all t > 0, (3.7) 

with 

9 — 72—1/ i\l O- 

C n+1 = C n {{C N /{2a)) l ' a ) 2 [T 2 - n -^ 
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This is inequality (13.11) for any p = 2 n+1 with n G N. 

We leave to the reader the proof that limsup n ^ 00 C n < oo. Hence, passing to the limit 
n — > oo in (13.71) we obtain inequality (13.11) for p = oo. 

The Holder inequality 

II || ^ || ||2 n+1 /p-l|| ||2-2 n+1 /p 
||^ Hp _ ||^||2 n Il^ll2 n + 1 

completes the proof for every p G (2 n , 2 n+1 ). □ 

Lemma 3.2. Let a G (1,2). Assume thatw = w(x,t) is the smooth approximation of the 
rarefaction wave, namely, the solution of problem {\2.2U^ - {\2.21^ . Then for each t > we 
have 

/■oo 

/ ||waxe(0llp dt < oo for every p G (1, oo] 
J to 

and 

[ t \\A a w(t)\\ p dt<Clog(2 + t) for p = ^, 

J to oc - 1 

all t >t and C > independent oft. 

Proof. It follows from the decay estimates recalled in Lemma 12.41 that 

/•OO poo 

/ \\w xx (t)\\ p dt<C r 3/a+i/(2p) rft < oo for every p G ( 1; oo]. 

By the interpolation inequality ( 12. 5ft and Lemma 12.41 we obtain 

||A Q «;(t)|| p < c(l+t)(- 1+1 ^( 2 - a )(l + t) M/2+1/(2p))(a_1) 
= C(l +t)-( 1+Q )/ 2 +( 3 - Q )/( 2 p). 

Hence, the rate of decay on the right hand side is equal to —1 for p = (3 — a) /(a — 1). □ 

Lemma 3.3. Let a G (1,2). Assume that u = u(x,t) is the solution of (I T. 1.2\) and 
w = w(x,t) - of (jOa> - (jQl . Suppose that u - w G L P (R) for p = (3 - a)/(a - 1). 

\\u(t) - w(t)\\ p < CTog(2 + t). 
Proof. Denoting v = u — w, we see that this new function satisfies 

v t + A a v + ^[v 2 + 2vw} x = -A a w + w xx . 

We multiply this equation by |v| p_2 v and we integrate over R to obtain 

\jt j \ v \ Pdx + J ( A "^)(l^r _2 ^) + ^ J [v 2 + 2vw} x \v\ p ~ 2 vdx 
(-A a w + w xx )(\v\ p - 2 v)dx. 
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(3i 



It follows from Lemma [2.31 that j R (A a v)(\v \ p ~ 2 v )dx > 0. The third term on the left hand 
side of (13. 8p is nonnegative by Lemma 12.41 and the same argument as the one used in the 
proof of Lemma I37T} cf. identity (13. 3p . Moreover, using the Holder inequality, we have 



(-A a w + w xx )(\v\ p - 2 v) dx 



< (\\A a w\\ p + \\w x 



ip-i 



xx\\p) 



Consequently, (I3.8P implies the following differential inequality 



|lK*)ll£<p(l|A a ™WiiP 



+ \\w xx 



XX \" ) 1 1 p 



p-1 

p > 



which, after integration, leads to 



K*)IIp< IK*o)H P + 



|A Q w(r)|L + \\w xx (r)\LdT. 



to 



The proof is completed by the result stated in Lemma 



□ 



Now, we are in a position to prove the main result of this paper 

The proof of Theorem First, we consider the auxiliary solution u = u(x,t) of the 
fractal Burgers equation (11 .11) with the step-like initial condition (11.101) . In this case, the 
measure m = (u + — U-)5q is nonnegative, hence by Theorem 12.21 u x > and by Lemma 

\\u{t) - u(t)\\ p < cr^/^lK - wolli 

for every p G [1, oo] and all t > 0. 

Next, we compare u with the smooth approximation of the rarefaction wave that is 
with the solution w = w(x,t) of (l2.20l) -f f2T2TT) (note that uo = wo). By Theorem 12.21 and 
Lemma [2 A\ we obtain 

II5 3S (0IL + Htw^lU < cr 1 . 



Moreover, using the following Gagliardo-Nirenberg inequality 

Ml^C^UM 1 -^ (3.9) 

valid for any 1 < p < p < oo and a = (l/po — l/p)/(l + l/po), we have 

\\u(t)-w(t)\\ p < C7(||5 B (t)|| 00 + \\w x (t)\U a \Wt) - w(t)\\l; a 
<Ct- a \\u(t)-w(t)\\l; a . 

Choosing p = (3 — a) /(a — 1) (hence a = [a — 1 — (3 — a)/p]/2), by Lemma 1331 we 
conclude \\u(t) — w(t)\\ p < Ct~ a log(2 + t) for every p G (po, oo]. Here, we are allowed to 
use Lemma [3731 because u — w G L 1 (IR) n L°°(M) C L P (R) for every p G [1, oo]. 

Finally, it follows from Lemma [2741 that the large time asymptotics of w(t) is described 
in L P (IR) by the rarefaction wave w R (x,t). 

The proof is complete because for 1 < a < 2 we have (1 — l/p)/a > (1 — l/p)/2. 
Moreover, since 1 < p < p, we have (1 — l/p)/2 > (l/po — l/p)/(l + l/po)- □ 
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4 Additional comments and possible generalizations 



Our main result is stated and shown in the simplest case of equation (11. II) . however, 
several generalizations are possible. 

First of all, the operator A a can be replaced by the Levy operator £ which is a 
pseudodifferential operator defined by the symbol a = a(£) > 0, Cv (£) = a(f)u(f). Here, 
the function e~ ta ^ should be positive-definite, so the symbol a(£) can be represented by 
the Levy-Khint chine formula in the Fourier variables 

a(£) = ibi + Q(0 + jf (l - e-** - itf * { \ v \<i}{v)) Wv)- (4.1) 

Here, b G K. is fixed, <3(£) — q£, 2 with some q > 0, and n is a Borel measure such that 
n({0}) = and J R min(l, \r]\ 2 ) U(dr]) < oo. 

Detailed analysis of conservation laws with the anomalous diffusion operator C is 
contained in papers [31 HI [5]. Here, we would like to emphasize that the fundamental nature 
of the operator £ is clear from the perspective of probability theory. It represents the most 
general form of generator of a stochastically continuous Markov process with independent 
and stationary increments. This fact was our basic motivation for the development of the 
theory presented above. 

In order to show the convergence toward rarefaction waves of solutions of conservation 
laws with the Levy operator, we need the counterparts of estimates (I2.13I) - (I2.14I) of the 
semigroup of linear operators e~ tc generated by — C They are valid e.g. under the 
assumption that the symbol a of £ has the form 

a(0=m a + H0, (4-2) 
where £>0,0<a<2 and A; is a symbol of another Levy operator K such that 

lim ^ = 0. (4.3) 

The assumptions f!4.2[) and (14. 3 p are fulfilled, for example, by multifractal diffusion oper- 
ators 

k 

C = -a dl + ^A-dl)^ /2 

with a > 0, aj > 0, 1 < ctj < 2, and a = mmi<j< k aj. We refer the reader to [HE] for the 
reasoning leading to the decay estimates of solution of nonlinear problem with operator 
L satisfying (I4.2l) - (l4.3p . That argument can be directly adapted to obtain counterparts of 
Theorem 12.21 and Lemma \3. II with A a replaced by C. Note here that the L p -L q estimates 
of the semigroup e~ tc are equivalent to a certain Nash inequality, see the papers [T5], H] 
and the references therein. 

Our result also holds true, if we replace the nonlinear term uu x in (11.11) by f{u) x 
with a strictly convex C 2 -function / (as in the paper of Il'in and Oleinik [10]) satisfying 
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f"(u) > k for some fixed k > and all u G KL Under this assumption, we immediately 
generalize Theorem 12.21 and we obtain the decay estimate (12. 3p . 

In order to show the counterpart of Lemma |3"TT| we should use the assumption f"(u) > 
k and to replace equalities (13.31) by the following (recall that v = u — u) 

[ [f(u) - f(u)] x \v\ P ~ 2 v dx > « (l - / u x \v\ P dx > 0. 

This argument, however, is known and used systematically e.g. in [22j inequality (3.5)] 
(see, also [101 [13 UHl EI]), hence we skip other details. 
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